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Abstract.  Using  observations  of  both  hard  X-rays  and  y-rays  in  the 
large  solar  flare  on  June  7,  1980,  we  infer  the  amount  of  chromospheric 
heating  due  to  bombardment  both  by  non-thermal  electrons  and  by  protons, 
respectively.  If  a  thick-target  model  for  the  X-ray  bremsstrahlung  is 
adopted,  then  proton  heating  is  shown  to  be  important  only  in  the  lower 
chromosphere;  however,  if  the  hard  X-rays  are  substantially  thermal  in 
origin,  then  proton  heating  may  play  an  important  or  indeed  dominant 
role  in  determining  the  structure  of  the  entire  flaring  chromosphere. 

1.  Introduction 

Observations  of  solar  flares  in  hard  X-rays  and  y-rays  provide 
information  on  the  populations  of  energetic  electrons  and  protons, 
respectively,  in  the  flare.  In  the  case  of  electrons,  the  observed  hard 
X-ray  emission  is  simply  proportional  to  the  rate  of  injection  of  energetic 
electrons  (except  for  very  rapid  events — Emslie,  1982),  and  one  can,  from 
pulse  height  bremsstrahlung  spectra,  deduce  the  form  of  the  injected 
energetic  electron  energy  distribution  at  each  instant  throughout  the 
flare.  The  determination  of  the  instantaneous  injected  proton  spectrum 
is  not  so  straightforward,  since  y-ray  emission  takes  the  form  of  line 
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emission  superimposed  cn  a  continuum,  and  since  some  y-ray  emissions 
(such  as  the  deuterium  formation  2.223  MeV  line:  Ramaty,  1982)  are  delayed 
phenomena.  Nevertheless,  it  is  still  possible  to  obtain  an  estimate  for 
the  time- integrated  flux,  or  fluence,  of  both  electrons  and  protons  from 
observation. 

Once  the  form  of  the  fluence  energy  spectrum  is  knovm,  it  is  a 
relatively  straightforward  matter  to  compute  the  variation  of  deposited 
energy  with  depth  throughout  the  atmosphere,  caused  by  the  collisional 
degradation  of  the  energetic  particles  (Emslie,  1978).  It  is  therefore 
possible,  and  of  course  of  interest,  to  compare  the  heating  rates  due 
to  both  electron  and  proton  bombardment  in  events  where  both  hard  X-ray 
and  y-ray  observations  are  available.  Our  conclusions  will  necessarily 
depend  on  Interaction  models  adopted  for  the  hard  X-ray  and  y-ray 
emission  processes,  since  this  affects  the  emitted  photon/inj ected  particle 
ratio.  Ramaty  (1982),  on  the  basis  of  a  comparison  of  interplanetary 
proton  spectra  and  y-ray  fluxes  in  a  number  of  events,  has  shown  that 
y-rays  are  almost  certainly  formed  by  a  thick  process  involving  non- 
thermal  proton  streams.  On  the  other  hand,  there  are  at  present  two  viable 
possibilities  for  the  hard  X-ray  emission  process — either  thick  target 
(Brown,  1971)  or  thermal  (Brown  et  al.,  1979;  Smith  and  Lilliequist,  1979; 
Smith  and  Auer,  1980).  Both  models  appear  to  be  consistent  with  hard  X-ray 
intensities  and  spectra  (Brown,  1974;  Elcan,  1978),  and  with  available 
polarization  measurements  (Emslie  and  Brown,  1980).  However,  a  thermal 
interpretation  of  the  hard  X-ray  burst  implies  fewer  electrons  precipitat¬ 
ing  into  the  lower  atmosphere  (Emslie  and  Vlahos,  1980;  Smith  and  Brown, 
1980),  and  so  a  lower  chromospheric  heating  contribution  from  electrons. 

In  this  article  we  therefore  compare  electron  and  proton  collisional  heating 
rates  under  both  interpretations  of  the  hard  X-ray  emission.  We  apply  our 


results  to  the  event  of  June  7,  1980,  for  which  good  observations  in  both 
hard  X-rays  and  y-rays  were  obtained. 
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2.  Observations 


At  approximately  0312  UT  on  June  7,  1980,  an  intense  solar  flare 
(X-ray  class  M7)  was  observed  near  the  west  limb  of  the  solar  disk  (coordin¬ 
ates  N12°  W74°)  in  active  region  2495.  The  event  is  characterized  by  a 
series  of  sharp  spikes,  observed  in  hard  X-rays,  y-rays,  and  microwaves, 
and  presents  many  fascinating  problems  of  interpretation  (Bai  et  al.,  1982; 
Kane  et  al. ,  1982;  Kiplinger  et^  al . ,  1982).  As  mentioned  in  Section  1, 
determination  of  the  instantaneous  injected  fluxes  of  both  electrons  and 
protons  is  not  possible.  However,  under  an  assumed  Interaction  model,  the 
particle  number  spectra  may  be  obtained  for  both  sets  of  particles,  and  this 
is  sufficient  to  compare  the  specific  energy  deposited  (erg  per  cm  of  vertical 
distance  per  unit  ambient  density)  by  electrons  and  protons.  In  Table  I  we 
present  the  parameters  of  the  particle  number  spectra  dN/dE  for  the  June  7, 
1980  event,  under  the  assumption  of  thick-target  models  for  both  hard  X-ray 
and  y-ray  production,  where  the  functional  form 


dN(E) 

dE 


(6-1) 


(1) 


has  been  assumed.  Here  N(E^)  is  the  total  number  of  particles  injected  with 

* 

energies  greater  than  E^  (arbitrary,  except  for  the  condition  E^  >  E  ;  here 

we  take  Ej^  ■  20  keV  for  electrons  and  10  MeV  for  protons) ,  6  is  a  constant 

* 

spectral  index,  and  E  is  a  lower  cutoff  energy  to  the  particle  number 
spectrum.  Note  that  the  spectral  form  (1)  is  principally  for  mathematical 
convenience  (see  Section  3),  and  should  not  be  taken  as  accurately  corres¬ 
ponding  to  the  actual  particle  spectra.  Also  note  that  the  time  and  area 
integrations  performed  in  arriving  at  the  values  in  Table  I  of  course 
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obscure  many  of  the  interesting  features  of  the  June  7  event  (e.g.  the 
multiple  spike  character,  also  possibly  related  to  a  multiple-loop 
character  for  the  event  {Emslie,  1981a]);  however,  this  will  not  affect 
the  order -of -magniturde  conclusions  presented  here. 

3.  Variation  of  Deposited  Energy  with  Depth 

The  transport  and  collisional  degradation  of  energetic  charged 
particles  through  the  solar  atmosphere  diring  flares  has  been  studied 
by  a  number  of  authors,  and  energy  deposition  rates  as  a  function  of 
depth  evaluated  (see,  e.g.,  Najita  and  Orrall,  1970;  Brown,  1972,  1973; 

Orrall  and  Zirker,  1976;  Emslie,  1978,  1980;  Hudson  and  Dwivedi,  1982). 

While  it  is  generally  accepted  that  energetic  electrons  may  be  responsible 
for  depositing  a  substantial  amount  of  energy  in  the  chromosphere  during 
the  flare  flash  phase,  the  role  of  protons  in  heating  the  lower  atmosphere 
is  still  somewhat  ct, 'troversial.  Most  of  the  discussion  regarding  proton 
heating  has  concentrated  on  the  lower  chromosphere/ temperature  minimum/ 
upper  photosphere  regions  of  the  atmosphere,  due  to  the  greater  penetrat¬ 
ing  power  of  deka-MeV  protons  compared  with  deka-keV  electrons  (Machado 
et  al. ,  1978;  Brown  et  al. ,  1981;  Hudson  and  Dwivedi,  1982),  but  few 
authors  have  considered  the  role  of  proton  heating  in  the  upper  chromo¬ 
sphere  (see  Lin  and  Hudson,  1976).  In  this  section  we  shall  follow  the 
analysis  of  Emslie  (1978)  in  order  to  assess  the  relative  roles  of  proton 
and  electron  heating  at  various  levels  of  the  atmosphere  in  the  June  7,  1980 
event . 

-3  -1 

The  energy  deposition  (erg  cm  s  )  due  to  the  passage  of  a  beam  of 

c 

charged  particles  through  a  relatively  cold  background  target  is  given 
by  Equation  (36)  of  Emslie  (1978);  using  the  spectrum  (1)  and  integrating 
this  rate  over  the  area  of  the  flare  and  over  the  duration  of  the  event  gives 
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the  total  energy  deposited  per  cm  of  height  throughout  the  event: 


where  K  ■  (e  *  electronic  charge  [e.s.u.J),  n  (cm  is  the  ambient 

density,  iiQ  is  the  cosine  of  the  injected  pitch  angle  (to  the  downward 

_2 

pointing  magnetic  field)  of  the  beam,  N(cm  )  is  the  overlying  proton  column 
density  at  the  point  in  question,  B  is  the  beta  function  and  the  para¬ 
meters  S  and  y  are  defined  as  follows  (Emslie,  1981b): 


.  2  *  A  + _(l-x)  A"  g 
A'  +  x  (A-A*)  ’  PP 


0 


m 

Ye  -  [x  A  +  (1-x)  A']  ,  Yp  *  ^  [x  A  +  (1-x)  A']  .  (3) 

e 


In  these  expressions  subscripts  e  and  p  refer  to  electron  and  proton 
bombardment  respectively,  x  is  the  fractional  ionization  level  of  the 


target,  mfi  and  nip  are  the  electron  and  proton  masses,  and  A,  A'  and  A" 
are  Coulomb  logarithms,  defined  in  Emslie  (1978).  (Note  that  the 


expressions  (3)  correct  an  error  in  Emslie  Cl 978 D »  see  Emslie  [l981b  J) 


Equation  (2)  is  valid  only  if  we  are  considering  values  of  N  > 

^E*2/ (2  +  -|)  ykJ,  where  E  is  the  cutoff  energy  in  the  injected  particle 
spectrum  (see  Equation  (1)).  For  values  of  N  less  than  this,  the  beta 
function  in  Equation  (2)  should  be  replaced  by  an  incomplete  beta  function 
(see  Equation  (8)  below),  which  results  in  a  reduced  heating  rate  Ig.  In 

*  V  * 

view  of  &he  uncertainty  in  the  value  of  E  ,  we  (formally)  take  E  as  zero 

20  -2 

for  both  electrons  and  protons;  for  N  =  10  cm  (see  discussion  in 

*  < 

Section  4),  Equation  (2)  then  yields  correct  values  for  Ifi  if  E^  ^  30  keV 
and  E  £  2  MeV,  which  seem  to  be  reasonable  conditions  to  satisfy 
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FIGURE  1 


(e.g.  Kane  et  al.  1980;  Ramaty,  1982). 

Since  we  are  here  interested  principally  in  the  ratio  of  energies 
deposited  by  proton  and  electron  bombardment,  we  can  take  the  target  to 
be  fully  ionized;  setting  x  -  1  in  equations  (3)  and  substituting  in 
Equation  (2)  results  in  the  required  energy  deposition  amounts: 


where  we  have  set  ■  1  (vertical  injection).  If  we  now  set  e* 

20  keV,  EijP  -  10  MeV,  and  use  the  values  of  and  6  from  Table  I,  we 
arrive  at  the  energy  deposition  amount  versus  column  depth  curves  of 
Figure  1,  to  be  discussed  in  the  next  Section. 

4.  Discussion 

From  Figure  1  it  Is  Immediately  apparent  that  electron  heating  pre¬ 
dominates  down  to  N  =  2.5  x  1022  cm-2  in  the  June  7,  1980  event,  at 
least  in  a  thick  target  interpretation  for  the  hard  X-ray  burst.  The 
level  of  the  flaring  chromosphere  lies  around  N  -  102®  -  102*  cm-2 
(Machado  and  Linsky,  1975;  Machado  et  al,,  1980);  thus  we  see  that  in 
such  an  interpretation,  proton  heating  is  unimportant  in  the  upper 
chromospheric  energy  balance.  We  note  in  passing  that  the  June  7,  1980 
event  has  a  remarkably  high  y-ray  to  hard  X-ray  flux  ratio 
(A.  L.  Kiplinger,  private  communication),  so  that  in  other  events  proton 
heating  will  be  of  even  lesser  importance. 

However,  in  a  thermal  interpretation  of  the  hard  X-ray  emission,  the 
situation  could  change  substantially.  As  remarked  in  Section  1,  such  an 
interpretation  involves  a  fewer  number  of  precipitating  electrons  (in 
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agreement  with  other  observational  determinations  of  this  precipitating 
flux  —  Brown  et  al.,  1978;  Emslie  et  al . ,  1978;  Donnelly  and  Kane,  1978), 
thus  reducing  the  value  of  Ig>e,  while  IgjP  remains  unaltered.  Vlahos  and 
Papadopoulos  (1979),  Emslie  and  Vlahos  (1980)  and  Smith  and  Brown  (1980) 
have  considered  the  precipitation  of  high  energy  electrons  from  a 
principally  thermal  source  at  a  temperature  of  several  x  10®  K,  with  some¬ 
what  varied  conclusions  regarding  the  role  of  such  precipitating  electrons 
in  the  flare  process.  Emslie  and  Vlahos  (1980)  consider  a  model  in  which 
the  beam  is  instantaneously  repopulated  by  direct  electron  field  accelera¬ 
tion  and  Dreicer  runaway  in  the  source;  the  injected  precipitating  flux  is 


F(E0) 


-E*/kT 

e 


electrons  cm 


-2  -1 

s  , 


E  2 
o 


(5) 


_3 

where  n  and  T  are  the  source  density  (cm  )  and  temperature  (K)  respectively, 
me  is  the  electron  mass,  k  is  Boltzmann's  constant,  vfi  *  (kT/mg)*5  is  the 
electron  thermal  velocity  and  E  is  the  minimum  energy  necessary  for 
escape  -  5kT  (Brown  et  al. ,  1979;  Smith  and  Brown,  1980).  Multiplying 
F(Eq)  by  AT,  where  A  is  the  injection  area  and  T  the  duration  of  the 
event  gives  N*,  the  number  of  injected  electrons  above  energy  E*. 

We  now  consider  the  following  source  parameters  (L  »  source  length) 


-  m11  -3 

n  -  10  cm 

T  *  3  x  108  K 

,  ,  .8 
L  -  10  cm 


A 


10 


15 


2 

cm 


T  -  60  s 


(6) 


6 


4 


2  45  -3 

appropriate  to  the  June  7  event  (with  emission  measure  n  V  *  10  cm  ) . 

* 

Substituting  these  values  in  eq.  (5)  gives  E  •  130  keV  and  Nj^q  keV  » 

3.6  x  1035;  eq.  (4)  then  gives  an  expression  for  the  energy  deposited  as 
a  function  of  depth  N  (Emslie,  1978;  Emslie  and  Vlahos,  1980): 


8  x  1049  N'2  a(N)  ;  N  £  1.7  x  1021  cm"2 
1.8  x  1050  N-2  ;  N  >  1.7  x  1021  cm“2. 


(7) 


where 

3KN 

*2 

a(N)  *7  x  (l-x)'2/3  dx-|  [l-(l+2>  (1-y)  /3],  (8) 

*2 

where  y  *  3KN/E  .  Note  that  the  energy  deposited  has  a  maximum  at 

y  =  1  (N  =  E  /3K),  since  no  electrons  of  energy  <  E  escape  from  the 

source.  For  small  N,  a  =  *j(3KN/E*2)2 ,  giving  I  /n  =  1.4  x  107  ergs, 

independent  of  N.  The  full  behavior  of  I„  /n  with  N  from  eq.  (7)  is 

21  -2 

shown  in  Figure  1.  We  see  that  below  N  -  10  cm  the  value  of  I 

“  > e 

(N)  is  not  substantially  different  from  its  value  in  the  thick  target 

case.  Physically  this  is  because  the  flux  of  high  energy  electrons  is 

similar  in  both  models.  However,  at  lower  N  values,  corresponding  to 

low  energy  electrons,  I_  in  the  thermal  model  is  much  smaller  than 

B,e 

in  the  non-thermal  case,  because  the  collisional  damping  rate  of  high 
energy  electrons  is  low,  and  there  are  no  low  energy  electrons  precip¬ 
itated  (as  these  would  be  in  the  thick  target  model).  Thus  for  N  £ 

20  -2 

5  x  10  cm  proton  heating  can  in  fact  dominate  over  electon  heating. 
This  value  of  N  is  entirely  consistent  with  the  position  of  the  upper 
chromosphere  in  flares,  especially  in  the  impulsive  phase  (Machado 
et  al. ,  1980;  Emslie  et  al . ,  1981),  and  so  we  conclude  that,  in  a 
thermal  interpretation  of  the  hard  X-ray  burst,  proton  heating  may  be 


important  in  the  chromospheric  energy  balance,  at  least  in  events  with 
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relatively  high  y-ray/hard  X-ray  flux  ratios,  such  as  the  June  7,  1980 
event.  Further  clarification  of  this  issue  must  await  a  better  under¬ 
standing  of  the  tail  repopulation  mechanism  in  hybrid  thermal/non-thermal 
models  (see  discussion  in  Smith  and  Brown,  1980) ,  and  an  anambiguous 
determination  of  the  thermal  or  non-thermal  character  of  the  hard  X-ray 
emission  in  flares. 
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Table  I 

Thick  Target  Injected  Number  Spectra 

* 

Electrons 

N,_  „  <S 

20  keV  e 

7  x  10  37  3.9 

*A.  L.  Kiplinger,  private  communication. 
tRamaty,  1982,  Table  3  and  Figure  1. 


Protonst 
N10  MeV 


32 

8.5  x  10  2.2 
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FIGURE  CAPTION 


Figure  1:  Energy  deposited  per  unit  ambient  particle  number  density  n 
by  electrons  (e)  and  protons  (p)  in  the  ^ne  7,  1980  event.  The  solid 
lines  refer  to  a  thick-target  model,  and  the  dashed  line  to  a  thermal/ 
thick-target  "hybrid"  model  of  Emslie  and  Vlahos  (1980);  see  Section  4. 


